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Models of spatially distributed populations and/or communities are still a 
matter of challenge for the students working in population biology, ecology, 
environmental sciences and mathematical modelling. An adequate model 
pattern to describe, model and predict the impact of spatial structure on 
a community dynamics, as well as the migration processes themselves, is 
a key problem here. Adequate modelling of a spatial transfer of a being 
is the basic difficulty here. Currently, the basic methodology addressing 
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the problem stands on implementation of partial differential equation of 
"reaction kinetics" type. 

Such approach has serious discrepancy. To be valid, the models strongly 
require that the beings move over space randomly and spontaneously (aim- 
lessly). This constraint is never met in nature; even microorganisms control 
their spatial redistribution [1-3]. 

Previously, there was proposed an approached to model a dynamics of a 
community with respect to spatial effects based on the (micro)evolutionary 
principle [3-5]. In brief, the principle forces beings to migrate in the manner 
improving their existence. An improvement of existence of beings is here 
the key question; the answer is given by the net reproduction function [5-8] 
k (p, r ). Here p = pi^r) is the (local) density of a population, and ~r is the 
point in space. Obviously, the population density p depends on the point at 
space. 

Net reproduction function results from two effects: the former is repro- 
duction, and the latter is inheritance [6-8]. If these two effects take place, 
then an equation of a dynamics of biological entity must look like 

p{-r 1 t)=p{~r 1 t)-k{p{~r,t)) or p t +i{~r , t) = p t ( r*, t) ■ k (p t ("r, t)) , (1) 

for discrete time. Here k(p) (k(N), respectively) is the net reproduction func- 
tion. It must be bounded above. The equations are the equations with 
inheritance. A comprehensive theory of such equations, in the most general 
case, including the investigation of three types of evolutionary stability see 
in [5-8]. The most general result is that net reproduction function k(p) must 
meet the extreme principle, for any free evolving biological community. k(p) 
is maximal for those species (entities) that sustained during the evolution; 
k(p) = in continuous time case, and k(p) = 1 for discrete time. 
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Finally, the basic issue of our model of spatially dependent dynamics of a 
community is that any migration must not be a random walk, but result in 
a growth of net reproduction. Further, we shall consider a model in discrete 
time and discrete space sites (called stations). A study of continuous model 
brings severe technical problems, so we shall start from a discrete case. 



1 Model of Two-Species Community 

We shall study a dynamics of a community consisting of two species; they 
are supposed to be "prey predator" related. It means, that one species 
exists due to external resources, but other one lives due to the beings of the 
fon„e r species. It is a.so supposed, that both spec.es occupy two station! 
and migration means a transfer of being (of any species) from station to 
station. Any other movements (inevitable in real situation) are neglected 
and supposed to have no effect on a community dynamics. 

The dynamics of (isolated) subcommunity occupying a station is sup- 
posed to follow the discrete analogue of classic Lotka-Volterra equation, if no 
migration takes place: 

N t+1 = N t -(a-bN t -fX t ) M t+1 = M t ■ (c - dM t - gY t ) 
X t+1 = X t ■ (efN t - hX t ) Y t+1 = Y t ■ (egM t - kY t ) . 

Here N t ,X t are the abundances of prey and predator, respectively, at the 
first station, and M t ,Y t are similar variables at the second station. Pa- 
rameters a and c determine a fertility of prey population, in the relevant 
stations; parameters b and d describe the density-dependent self-regulation 
of this population, in relevant station. Parameters h and k describe similar 
1 These are the sites considered together with the environmental conditions. 
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density-dependent regulation at predator subpopulations. Parameters / and 
g describe, in general form, an efficiency of the interaction of the beings of 
these two species, including a success of hunting, success of escape, etc., in 
corresponding stations. Finally, e represents an efficiency of the conversion 
of prey biomass into the predator biomass. 

1.1 Basic Model of Migration 

Parameter p, < p < 1 is a mobility of prey beings; similar, q, < q < 1 
is mobility of predator beings. These parameters are the transfer cost and 
might be interpreted as a probability of the successful migration from one 
station to other; success here means that no damage for further reproduction 
had taken place. Migration from station A to station B starts, if living 
conditions "there" are better, than "here", with respect to the transfer cost: 



for prey and predator beings, respectively. It should be stressed, that the mi- 
gration act is executed independently by each being, while the model consid- 
ers it as a population event. The backward migration conditions are defined 
similarly: 



Migration act runs each time moment t, for both species independently. If 
neither of the inequalities (pH HD are fulfilled, then no migration takes place, 
at the given time moment t. Prey migration flux A (predator migration flux 



(a-bN t -fX t ) < p .(c-dM t -gY t ) , 
(efN t - hX t ) < q ■ (sgM t - kY t ) , 



(3) 



p .(a-bN t -fX t ) > (c-dM t -gY t ) 
q ■ (efN t - hX t ) > (egM t - kY t ) . 



(4) 
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(5a) 



(5b) 



0, respectively) must equalize inequalities ([3111]): 

(a - b(N t - A) - fX t ) = P -(c- d(M t + pA) - gY t ) , 
(efN t - h(X t - 9)) = q ■ (egM t - k(Y t + qG)) 
for the case (j3J), or 

p • (o - 6(iV t + pA) - /X t ) = (c - d(M t - A) - ^FO , 
g • (e/JV t - h(X t + g9)) = (^M t - k(Y t - 6)) 

for the case (j3J). Then, A (0, respectively) is equal to 
^_pc-a + bN - pdM + fX - pgY _ hX + eqgM - efN - qkY 
b + p 2 d h + q 2 k 

(6a) 

for migration form station A to station B, and 
^_pa-c + dM - pbN + gY - pfX _ kY + eqfN - egM - qhX 
d + p 2 b k + q 2 h 

(6b) 

for the backward migration. 

Finally, let's outline how the basic model Ml works. For each time mo- 
ment t, a direction and the migration fluxes (A and G, respectively) are de- 
termined. Then, the species redistribute themselves according to the Eqs. flU}. 
Then, the abundances of the next generation {N t+ i, X t+ i, M t+ i,Y t+ i} are 
determined, according to (j2J), with the relevant abundances of the current 
generation {N t ,X t ; M t ,Y t } defined by (jHD- If no migration must take place 
at the current timer moment t, the the stage with species redistribution is 
omitted. 



1.2 Reflexive Behaviour 

Reflection in behaviour means an ability of a being to foresee and/or predict 
the behaviour of an opponent, in a competitive behavioural act. An imple- 
mentation of reflexive behavioural strategy by animals is a well known. Not 
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discussing here psychological or ethological aspects of such strategies imple- 
mentation, let concentrate on a simple model revealing the dynamic effects 
of them. 

Basic model (F2]-[6]) does not exhibit any reflexive behaviour. An intro- 
duction of that latter into the basic model may only be concerned with the 
spatial redistribution. In other words, reflection of the optimal migration be- 
haviour means that a being is able to "foresee" the migration behaviour of a 
competitive species being. With respect to it, one may assume the following 
patterns of the reflection in the behaviour of the species: (i) preys reflect 
predators; (ii) predators reflect preys, and, finally, (iii) both species reflect 
each other. 

Thus, within the framework of our model, a reflection means that the 
species manifesting a reflection in the behaviour, detects the migration con- 
ditions and chooses the migration flux according to the abundances of a 
competing species, that would be produced due to the migration of that lat- 
ter, not the current ones. In case (i) formula for would remain the same, 
but the formula for A would change for 



in dependence of the migration direction. Here X and Y are determined 
according to ([5]). 

Reciprocally, N and M are determined according to §5§, for the case (ii), 



A 



= < 



pc-a + bN -pdM + fX - pgY 
b + p 2 d 

pa-c + dM - pbN + gY - pfX 
d + p 2 b 



or 



(7) 
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but migration flux of predator would be determined by 

' hX + eqgM - efN - qkY 



e 



or 



h + g**_ (g) 

kY + eqfN - sgM - qhX 
k + q 2 h 

in dependance of the migration direction. 

Finally, if both species reciprocally reflect the behaviour of each other, 
then basic model should be changed for the following one. On the first stage, 
both species determine the migration fluxes according to basic model (0 [6]). 
Then, they redefine the migration fluxes (and migration direction, as well as 
the fact of migration) so, that each species changes the current abundances of 
the competitive beings for those that could be produced due to a migration 
rule determined by the basic model. So, they redefine the fluxes, redistribute 
themselves between the stations, and reproduce. 



2 Results and Discussion 

Main purpose of this paper is to figure out the sets of the parameters provid- 
ing an evolutionary advantage to a bearer of some (reflexive, or not) spatial 
distribution strategy. Evolutionary advantage here is understood as an excess 
of the total abundance of some species realizing reflexive strategy, in com- 
parison to the same species in case of realization of regular (non-reflexive) 
strategy. 

Table [I] shows the result of simulation observed for the following param- 
eters sets: 

#1 a = 3.1, c = 1.49999, b = 0.00098, d = 0.00099, h = 0.00052, k = 
0.0005, f = g = 0.00542, e = 0.099, p = q = 0.99. 
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#2 a = 2.5, c = 1.49999, b = 0.00098, d = 0.00099, h = k = 0.0005, 
f = g = 0.0054, e = 0.099, p = q = 0.99. 

#3 a = c = 1.5, & = d = 0.0001, /i = 0.0001, k = 0.0005, / = g = 0.0059, 
e = 0.05, p = q = 0.99. 

#4 a = c = 1.5, b = d = 0.0001, h = k = 0.00001, / = 0.0059, g = 0.0049, 
£ = 0.05, p = q = 0.99. 

#5 a = 1.7, c = 1.5, b = d = 0.0001, h = k = 0.00001, / = # = 0.005, 
e = 0.05, p = q = 0.99. 

#6 a = 2.2, c = 2.1, b = d = 0.0001, h = k = 0.00001, / = # = 0.005, 
e = 0.05, p = q = 0.99. 

First of all, it should be said, that the basic model ([2] - [H]) exhibits a 
great diversity of limit regimes. It may be a steady state (in both stations, 
for both species), limit cycles of various length, and a complex irregular 
behaviour looking like a dynamic chaos. All these peculiar regimes may be 
met in combinations, with respect to a station and/of a species. In general, a 
decrease of transfer cost p and q yields a simplification of an observed regime. 
Both the basic model (T2] - [6]) and its versions implementing various reflexive 
strategies of spatial distribution yield an expansion of the area of permissible 
parameter values, and the area of the phase space (i.e., abundance figures). 
All these issues are very interesting, from the point of view of the study of the 
models of optimally migrating communities, but they fall beyond the scope 
of our research. 

What we do, was a comparative study of those four models (basic model 
and three versions with reflexive behaviour) from the point of view of the 
evolutionary advantage. In fact, we tried various combinations of the pa- 
rameters, similar for all four models, in order to identify the model that 
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Table 1: Comparison of various strategies of space distribution. S — type 
of strategy: 1 - basic model, 2 - model i, 3 - model ii, 4 - model Hi. 



s 


N 


X 


M 


Y 


S 


iV 


X 


M 


Y 


Set of paramenter # [TJ 


Set of paramenter # [2] 


1 


2089 


5 


4539 


6 


1 


1513.5 


1.4 


1952.5 


56.2 


2 


1930 


20.5 


2368 


50 


2 


1438.8 


4.1 


2244.6 


55.6 


3 


1692.9 


43.7 


1849.9 


41.7 


3 


1524.5 


0.47 


1938.6 


56.7 


4 


2142 


0.1 


5040 


0.1 


4 


1530 


0.01 


5047 


0.01 


Set of paramenter # [3] 


Set of paramenter # [4] 


1 


3697 


18 


3698 


18 


1 


3542 


20.9 


4073 


14.3 


2 


4999 


0.01 


4999 


0.01 


2 


4997 


0.01 


4998 


0.01 


3 


3873 


6.4 


3869 


6.5 


3 


3619 


11.2 


3739 


8 


4 


4189 


13.8 


4190 


13.8 


4 


4325 


12.5 


4384 


11.5 


Set of paramenter # [5] 


Set of paramenter # [6] 


1 


6994 


0.05 


4996 


0.01 


1 


5.4 


0.1 


4.4 


0.04 


2 


5688 


32.1 


4122 


11.6 


2 


8406 


80.5 


7172 


69.3 


3 


4664 


19.5 


2961 


13.9 


3 


3.6 


0.07 


3.6 


0.1 


4 


5412 


30.5 


4086 


18.7 


4 


8404 


80.5 


7171 


69.3 



yields the highest total abundance of the species. 

To answer this question, we have carried out a series of simulations 
(computational experiments). We calculated the abundance of each species, 
in each station, for four models with the same parameter set. Then, the 
abundances of prey subpopulation (of predator subpopulation, in turn) were 
added. Table [1] shows the results of such comparison. The parameters yield- 
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ing the observed regimes are shown below the Table. 

It is evident, that there exist the parameters sets yielding an evolutionary 
advantage for various types of space distribution strategies. Yet, we did 
not studied carefully the peculiarities of the limit regimes relevant to each 
evolutionary advantageous situation, meanwhile, one may expect that the 
reflexive strategies provide an advantage for rather regular limit regimes, 
while the non-reflexive strategy of space distribution is advantageous for 
chaotic-like, complicated limit regimes. 
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